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Abstract 



We analyze the dynamics of a qubit-resonator system coupled with 
a thermal bath and external electromagnetic fields. Using the evolu- 
tion equations for the set of Heisenberg operators, that describe the 
whole system, we derive an expression for the resonator field, account- 
ing for the resonator-drive,-bath, and -qubit interaction. The renor- 
malization of the resonator frequency, caused by the qubit-resonator 
interaction, is accounted for. Using solutions for the resonator field, 
we derive the equation describing qubit dynamics. The influence of 
the qubit evolution during the measurement time on the fidelity of 
a single-shot measurement is studied. The relation between the fi- 
delity and measurement time is shown explicitly. Also, an expression 
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describing relaxation of the superposition qubit state towards its sta- 
tionary value is derived. The possibility of controlling this state, by 
varying the amplitude and frequency of drive, is shown. 

1 Introduction 

The possibility of achieving quantum coherence in macroscopic Josephson 
junction (JJ) circuits, envisioned by Leggett in the early 1980s P]-[3], was 
demonstrated experimentally almost 20 years later by Nakamura et al. [I]. 
Now, superconductive qubits are promising building blocks for the realization 
of a quantum computer. Progress towards quantum computing depends on 
the development of measurement schemes and readout devices. Each read- 
out process requires a finite-time interval during which the system evolves 
under the effect of the measuring device [S]-[7] and uncontrollable external 
j8],[9] or intrinsic [10]- [H] noises. This requires improving the quality of the 
Josephson contacts and effective isolation of qubit systems from the "elec- 
tromagnetic environment". There is stable progress toward the solution of 
the first problems. (See, for example, Ref. [15J.) The second problem is 
fundamental in nature because, for example, it is impossible to isolate the 
qubit system from vacuum zero-point oscillations or to avoid the influence of 
the measurement device on qubit relaxation or decoherence. 

In this paper, we study the simultaneous influence of external electro- 
magnetic fields and a thermostat (bath) on the qubit. It is assumed that 
this influence is not direct but via a resonator which is weakly coupled with 
the qubit. Therefore, the resonator response on the microwave field depends 
on the qubit state. In particular, by measuring the phase of a microwave 
field reflected from the resonator, one can perform a nondemolition mea- 
surement of the qubit state. The theoretical description of this scheme and 
its practical implementation are elucidated in numerous publications [2]- [I], 
[TB]-[20j. At the same time, it should be emphasized that this kind of mea- 
surement results in both the dephasing of the qubit wave function and partial 
relaxation caused by the effects of drive and the thermostat during the mea- 
surement. These effects decrease the measurement fidelity. By studying the 
qubit dynamics, we are able to estimate the influence of drive and bath on 
the fidelity. 

Our further consideration is based on the equations of motion of the 
Heisenberg operators. These operators describe the dynamics of both the 
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qubit and resonator field. The solution of these dynamical equations will 
allow us to estimate the upper limit of the measurement fidelity. 



2 Model 

We consider the system described by the following Hamiltonian: 
H = --tku q a z + %u r [a + a + + J2fhu n (b+b n + - 

igo y {a+ - a)j +ih^f n (b n a + - b+a^J + ihf (ca + - c + aj , (1) 

where the first, second, and third terms in the right side are Hamiltonians 
of the noninteracting qubit, the resonator, and the bath, respectively. The 
quantity, cr z , is the Pauli operator. We use here the diagonal representation 
for Hamiltonian of the isolated qubit. Different types of superconducting 
qubits, for example, charge qubit [5]- [7] or flux-biased phase qubit [13] can 
be used in practice. 

The quantities a + , 6+, c + are the creation operators of the corresponding 
excitations and a, b n , c are the annihilation operators. Usually, the drive 
variables, c + and c, are considered to be the classical quantities, c* and c, 
with given dependences on time: c*, c ~ e ±Wdt in which uj^ is the frequency 
of microwave field. 

The fourth, fifth, and sixth terms are the resonator-qubit, -bath, and 
-drive ineractions, respectively. The constants g, f n , and fo, describe the 
corresponding interaction strengths. u q is the transition frequency between 
qubit levels. The quantities, u r and ui n , are the frequencies of the resonator 
(cavity) and bath oscillators, respectively. The resonator Hamiltonian de- 
scribes the parallel connection of the local capacitance, C, and the induc- 
tance, L (co r = (LC)- 1/2 ). At the same time this lumped-element circuit 
can represent a transmission (microstrip) line if we deal with signals whose 
characteristic frequencies are sufficiently close to one of the eigenfrequencies, 
u r , of the line. More details can be obtained, for example, in [S] and [22] . 

Similar to papers [21], [22], we model the bath as an infinite set of har- 
monic oscillators with frequencies, u n . It can be seen from the explicit 
resonator-bath interaction Hamiltonian that the presence of the bath is im- 
portant even at zero temperature. In this case, only terms describing the 
annihilation of resonator excitations and the creation of bath excitations 
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appear. At finite bath temperatures, there are fluxes of enegy in both direc- 
tions: into and out of the bath. There is no unique description of the effect of 
the bath on the qubit-resonator system. For example, the model of a direct 
qubit-bath interaction is used in Ref. [23J. 

The resonator-bath as well as the resonator-drive interactions affect the 
qubit state due to the qubit-resonator interaction. This situation is very 
similar to the case considered in [5] in which semiclassical noise of the biased 
current causes qubit decoherence and relaxation. 



3 Evolution of the resonator field 

Using the Heisenberg representation for the operators in Eq. [U we can 
express the time derivative of a as 

a = -^{a,H] = -iu r a + jO y + fnK + foe (2) 

Here, the dependences of a y and b n on time are not yet specified. To de- 
termine the explicit dependences of u y (t) and b n (t), we will use an iterative 
procedure that assumes that all interaction parameters are small quantities. 
As in Eq. [2] we have 

d y = -u q a x , (3) 

I' 

From these equations we can obtain relationships for the more convenient 
variables a± = \(a x ± i^y)'- 

a± = ^iujq<j± + i^a z (a + - a). (4) 

Eqs. H] can be rewritten in the equivalent form as 

rt 



d x = LOqUy + i2—a z (a + — a). 



a±(t) = a ± (t )e^"^ to) + d [ dt'e T ^ (t - t ' ) [a z (a + - a)] t ,. (5) 



Further analysis is facilitated by assuming the qubit-resonator detuning is 
small: 



|^gr| ^qi ^ri ^qr — ^q — 0J r . 
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(We assume this inequality is satisfied throughout the paper.) At the same 
time the detuning should be large compared to the interaction parameter g : 

«1. (6) 



H | OJqr | 



In many papers (see, for example, Refs. [5]- [7]), the quantity in the left side 
of Eq. [6] is also considered as small parameter. In this case, the dynamics of 
the system can be analysed using a simpler (renormalized) Hamiltonian. 

Considering a z (t')a + (t')e~ lbJrt and <r z (t')a(t')e tulrt as slowly varying func- 
tions of t', we can integrate over t' in Eq. [5j The result is 

ll 0Jq r 

a_(t) = o_{t) - , 

where a± = a±(to)e TluJq( - t ~ to \ In the course of integration, we have assumed 
that t — to — > oo. 

Using Eqs. [7] we can express cr y (t), which appears in Eq. [21 in terms of 
a±(t),a z (t),a + (t), and a(t). The dependences of b^,b n on time can be also 
expressed in terms of a + and a. As in Eqs. HI we have 

b n = -iu n b n - f n a, (8) 

or 

b n {t) = b n {t) - f n [* dte^Watf), (9) 

J to 

where ~b n {t) = b n (t )e- iuln{t - to) . 

Multiplying both sides of Eq. [9] by f n and summing over n we obtain 

Y,fnh(t) = Y,fnbn(t) ~ f dt £ ^V^"'^') • (10) 
n n *0 n 

Using a simple approximation for the sum in the integrand of Eq. [10] 

n 

(see Refs. [2T],[22j), we can rewrite Eq. [10] in a simple form: 

E/»M*) = £M»(*)-£°(*)- ( 12 ) 



Then using Eqs. [7]and[T2J the equation for a(t) reduces to 



a = f c + fnb n - «| ( cr+ - a-) - ix^ z a + , (13) 



where x = 9 2 /(h 2 u qr ). 

The quantity, u r — x <J z, is the resonator frequency renormalized by the 
qubit-resonator interaction. This renormalization effect can be derived from 
the Jaynes-Cummings Hamiltonian using a unitary transformation, assuming 
that g/(hoj qr ) is a small parameter. 

The parameter, k/2, appearing in the bath- resonator interaction, is an 
important characteristic of the resonator. It can be seen from the structure 
of Eq. [13] that k/2 describes the field dissipation caused by this interaction. 
The ratio u r /n is the resonator quality factor, Q. Usually, high-Q resonators 
are used for qubit measurements. 

The solution of Eq. [HI] (in which the influence of the initial condition or 
transient stage is ignored) is given by: 

a(t) = iMt \ + V lfX{t \ + f " +(t) - , (14) 

LUdr + iK/2 „ UJ nr + %Kj2 nUgr + iK/2 

where u ir = uj-i — u r + x<r 2 , i = d,n, q. In the course of solution of Eq. [TBI 
explicit dependences of c, b n , a± on t were used. Also, the contribution of 
terms with cx_ and a + was neglected. This approximation is accurate when 

x, (g/h) « u r - 

The value for a + can be obtained from Eq. [lousing hermitian conjuga- 
tion. It is given by 

+/ t x _ ifpc + (t) ] J2 l f*fitM 9 cr-(t) (15) 

UJ r d + irl/2 ~ / UJ rn + iK,/2 U uj rq + m/2 

Using Eqs. dU [13 and [TTJ, we can easily show that the standard com- 
mutation relations between operators a and a + are fulfilled with an accuracy 
valid up to a small value of the order g 2 / (h oj 2 r ) if the drive variables are con- 
sidered as classical quantities. The deviation of [a, a + ] from unity is within 
the accuracy of perturbation procedure used here. 

It follows from Eqs. [TH and [TS] that the effect of bath is represented 
by the second ("noise") terms and the imaginary summand, in/2, in the 
denominators. The effect of drive on the resonator field critically depends 
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on the detuning Cj dr = uj d — u r + x°z- The field amplitude and the photon 
number in the cavity, n r = a + a, are the largest when the detuning is of the 
order of k/2. If the drive frequency is fixed, the detuning depends on the 
qubit state. For the upper and lower states, the resonant conditions can 
be very different when \ > K /2- This circumstance is commonly used for 
measuring qubit states by means of microwave fields (see, for example, Ref. 

inn- 

To estimate the importance of the different terms in Eqs. [14] and [151 
we will calculate the average (over bath variables) photon number in the 
resonator, n b r1 considering the qubit to be in the excited [a z = —1) or ground 
(<r z = 1) state. (The frequency, u dr , is equal to u d — u r =F X f° r °z — Tl, 
respectively.) Using Eq. [TT]and the relationship 

°-<r+ = ^ - a z), (16) 

we obtain 



n r ~ ~2 i , 2 I A (b n bn) ui n =ui r + ,2 ~2 ' (^) 



\foc\ 2 , v g 2 I-*. 



7,2 
qr 



where correlations between different bath modes were ignored ((&+&„/) ~ 
S ntn '). In the course of derivation of the second term in Eq. [T7|, we have 
considered that the average (b^b n ) depends on n via u n only. Then using 
Eq. [11] we were able to sum up over bath modes as: 

fn(bnbn) \ f°° j, , , ~, \ fn (^n ^n) \u) n =u+ui r -x<Tz 



UJr. 



(b+b )\ — r -V f1J{"-*nr)T = l h+h VI 

\ u n u n/\u) n =u) r / 2 1 2 //I / O ^ •> n \ u n u n/ \uj n =ui r - 

J— 00 uV T" rC / 4 J —00 Z7T 

The remarkable peculiarity of Eq. [17] is that the direct contribution of 
bath does not depend on the interaction constant, f n . This is in contrast to 
Eqs. [T41 and [151 for the resonator field. 

For thermal equilibrium, the average occupancy of the bath modes is 
given by the Bose-Einstain distribution function ube 
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(b+b n ) = n BE (u n ) = [e k BT _ i j (is) 

Therefore, in the absence of the drive and the qubit, the number of photons 
in the resonator is equal to that in the corresponding bath mode. In other 
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words, the bath and resonator temperatures are equal in this (equilibrium) 
case. 

Eq. [T7] shows explicitly when the driving field dominates the noisy in- 
fluence of the bath. Besides that, it follows from Eq. [T7] that the qubit 
"delivers" an almost negligible portion of photons to the resonator (the last 
term in Eq. H7|) even in the most favorable case, a z = —1. The physical 
reason for this is in the qubit-resonator detuning. The detuning decreases 
the probability of qubit excitations to "penetrate" into the cavity. 

Eqs. [14] and [15] can be used to study the fluctuations of photon numbers 
in the resonator. These fluctuations are responsible for qubit decoherence 
(more details can be obtained, for example, in Ref. [B]). 

In the next Section, we will use Eqs. [TJ] and [TS] to describe the qubit 
dynamics. Drive- and thermostat-induced variations of the qubit state during 
the measurement time will be studied. These variations are responsible for 
reducing the measurement fidelity of the qubit state. 

4 Qubit evolution 

The time variation of qubit states occupancies can be expressed in terms 
of the average value of the operator a z . When the qubit is in the state 
ip(t) = a(t)\0) + (3(t)\l), the average value of o z is: 

(a z }} t = (mw z \m} = Ht)\ 2 -\m\ 2 - 

If we know (<J z )t, the occupancies of the levels can be obtained from: 
l«W| 2 = ^(l + (^)t), |/3(t)| 2 = ^(l-K),). 

We will replace (& z )t by (cr z (t)), where a z (t) is defined in the Heisenberg 
representation and averaging is over the initial state, ip(to). Therefore, the 
evolution of the occupancies can be obtained from the temporal dependence 
of the operator a z (t). 

To obtain a z (t), we use, as previously, the equations of motion for the 
corresponding operators. The time variation for a z (t) is: 

& z = ^-[cr z ,ig(r y (a + - a)]. (19) 
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We can express a y (a + — a) as — i(a + a + + cr_a) in the spirit of rotating-wave 
approximation. Then Eq. [TjJ] reduces to 



■ 2 9, + ^ 
a z rti — 2— (a + cr — acr_). 



(20) 



In the next step, we will express the right-hand side of Eq. [20] in terms 
of a z {t). Similar to the consideration in the previous Section, we use the 
equation of motion for the operators a + a + , cr_a. Thus we have 



(d t + iuj qr )a + a + = i^a z a + (a + -a) a z ) + a + ( f c + + ^2 f n b+) . (21) 
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Using Eq. [12] and neglecting the term containing a + a + , we can rewrite Eq. 



as 
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d t + iuj qr + — )o"+a = —ij<r z \ a + a 



■*|-+^ + U° c+ + £^«)- (22) 



Using Eqs. [7J and HU we obtain from Eq. 



d t + iu qr + - ) cr + a + = a + ( f c + + J2 fJ>. 



■ 9 



n. \ a a + - 



+ 



9°z 



foe 



+ Ett 



fn"n 



f0C + + Y,fn>K> 



(23) 



fiuj qr \uJdr + in/2 „ uj nr + in/2, 

Let us average both sides of Eq. [23] over the thermostat variables. Then 
again considering slowly varying function of time, we can easily 

obtain (<7 + a + )b a th from Eq. [23] in the form 



(a+a + ) 



ig\e 



bath 



h(iui qr + K,/2) 



<t- ! n r + - 



U qr \UJdr + i^/2 „ u}^ + i«/2 



(24) 



where we have neglected the initial (at t — 0) correlations of the operator <t + 
with the operators a + ,c + . The value of is given by Eq. [17] in which the 
last term, representing the contribution of qubit, can be omitted due to its 
small value. 
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The expression for {aa z )b a th can be derived in a similar manner. Then 
the equation for a z is: 



<7z 



2g 2 
h 2 u 2 qr 



a 



2co qr sin(u qr t)e 2 * + /t( 1 — cos(co qr t)e 
\foc\ 2 



2u 



dr 



COS( 



(uj qr t)e 2 J 4- Ksin(u qr t)e 2 ' 



(25) 



Because u qr » k/2, we can omit in Eq. [25] the oscillating terms that are 
proportional to k and Cjdr- Thus, the rate equation reduces to 



2g 2 



h 2 u 2 qr 



2uj qr sin(uj qr t)e 2 +k 



2-^^ + 2 



a, 2 * M2 .y (26) 



In the absence of drive (/o = 0), Eq. [26] describes small- amplitude 
thermostat-induced Rabi oscillations with frequency u qr . These oscillations 
decay during the coherence time, ~ 2/k, of the resonator field. The oscilla- 
tions are accompanied by a slow qubit relaxation to the stationary value 



a st = (l+2(b+b n )\ Un -_ 



(27) 



(See, for example, Ref. [21] ■) The amplitude of oscillations depends on: 
the occupancies of the bath states, (&^& n ) U n =w r , the cavity losses, k/2, and 
the qubit-resonator detuning, u qr . In the case of an equilibrium bath, a st 
coincides with the equilibrium value known in the literature. It can be derived 
in an alternative manner using the density matrix formalism. 

The case of /o 7^ is more interesting in view of the possibility of using 
the drive to measure (control) the qubit. The preferred setup is realized if 
the microwave field is in resonance with the cavity-qubit system in which the 
qubit is in a given eigenstate (for example, in the excited state: Ud = u r + x) ■ 
(This is the case for which the last term in braces of Eq. [26] vanishes.) The 
measuring device can resolve the resonator frequencies w r ±x, corresponding 
to the different eigenstates of the qubit, only if the measurement time, r, 
is greater than (2%) _1 . At the same time, r should be as small as possi- 
ble to provide rapid control. Moreover, during the measurement time, the 
qubit evolves as described by the rate equation f [26|) . thus decreasing the 
measurement fidelity F. As a consequence of the measurement, an arbitrary 
superposition qubit state collapses to one of the eigenstates with a z = ±1. 
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The characteristic time of the collapse (defined here as the decoherence time) 
is assumed to be considerably shorter than r. Therefore, the fidelity of the 
measurement, F, for the qubit in the post-collapse state, |1 >, is defined as 

1-M 1/2 1 9 2 (n b r + l) 



T Jo T Jo 



dt\m = zi dt \—r) *i- ;v 7 (i 



qr 



(28) 

For simplicity, we have considered the case r < 2/k in which a z is given by 

— -sin 2 ' 

Taking into account Co> gr r >> 1, we have 

9 2 (n b r ±l) 
j 2 

qr 



1 - i^p. (29) 



This formula illustrates the effect of the Rabi oscillations, generated by the 
external drive and bath, on the fidelity. In contrast to cr z (t), the fidelity 
does not display the oscillating behavior. This is because the value of a z 
is averaged over the interval r >> |tu ?r | _1 . For longer measurement times, 
r > 2/k , all terms in Eq. [26] should be used to calculate a z (t) and F(r). 

It seems from Eq. [29] that the fidelity can be improved for smaller inter- 
action parameters. But decreasing g will decrease x = 9 2 1 (ft 2 ^qr)- m view 
of the inequality r > (2%) _1 and using Eq. [29] we obtain 

T > w^k- (30) 

It follows from Eq. [30] that as F approaches unity, the measurement time 
should be increased. 

We have described theoretically a single-shot measurement in which the 
initial qubit was assumed to be in either the excited or the ground state. The 
evolution of the superposed state can be investigated experimentally by the 
repetition of many single-shot measurements. The rate equation for (a z (t)) 
[the averaging is over the initial state, if)(t = 0)], corresponding to this kind 
of the experiment, can be derived in a manner similar to Eq. [26] It is given 

by 



tl UJ qr 



(a z )(cp + + <p-) + ip- -if + + 1 



(31) 
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where the initial (oscillatory) stage of evolution is ignored; 



^.. = ±1), vM = nl + l + ^JML_. (32) 
In the derivation of Eq. [31] the following identities were used: 
<PM = i[(l + a 2 )v9 + + (l-a 2 )^-]. 

and 

crzVivz) = ^K 1 + ff z)<P + ~ (1 - O^"]- ( 33 ) 
The solution of Eq. EH is given by 

M*)> = ^ + = 0)) - ^)e" 7 \ (34) 

in which the stationary value, a s t, and the relaxation constant, 7, are given 
by 

°st = — — I—, 7 = Z2 9 • ( 35 ) 

In the limiting case /o = 0, a st reduces to the previous result given by Eq. 
1271 In the limit of dominating drive in Eq. [17] (n r » 1 for the resonance 
conditions, Ud = u; r ± x), the qubit relaxes to the ground or excited states, 
respectively. In the case of large resonator-drive detuning, \ud r \ » X-> the 
last term in Eq. [32] can dominate. Then (a z ) relaxes to zero. 



5 Conclusion 

We have described a qubit-resonator system with an external drive and ther- 
mostat. Our consideration is based on the equations of motion of the op- 
erators in the Heisenberg representation. This is in contrast to the widely 
used density matrix approach. Considering the resonator-bath, resonator- 
drive, and resonator-qubit interactions as weak perturbations, we have de- 
rived expressions for the resonator field including the renormalization of the 
resonator frequency caused by the qubit-resonator interaction. A weak qubit- 
environment interaction is a necessary condition for reliable isolation of the 
qubit from the "external wold. 
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Also, we have derived the rate equation for the qubit variable, o~ z , describ- 
ing the occupancies of the qubit levels. The solution of this rate equation 
enables us to calculate the measurement fidelity, F, and to determine the 
dependence of F on the measurement time, r: increasing fidelity requires 
increasing the measurement time. (See Eq. [30j) Both quantities are very 
important parameters in view of the practical implementations of qubits. 
Therefore, the optimal choice of F and r should be carried out considering 
their interdependency given by inequality |30j 

The qubit relaxation, caused by the interactions with the bath and the 
drive, can be used to control the final qubit state. In particular, it follows 
from Eqs. I3"TH3"51 that, by varying the frequency and amplitude of the drive 
as well as the interaction time, t, we can get a qubit with a predetermined 
probabilities to be in the ground or excited state. 
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